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The large TV limit of SU(TV) gauge theories in 3+1 dimensions is investigated on the lattice by extrapolating 
results obtained for 2 < TV < 5. A numerical determination of the masses of the lowest-lying glueball states 
and of the topological susceptibility in the limit TV — > oo is provided. Ratios of the tensions of stable fc-strings 
over the tension of the fundamental string are investigated in various regimes and the results are compared with 
expectations based on several scenarios - in particular MQCD and Casimir scaling. While not conclusive at zero 
temperature in D=3+l, in the other cases investigated our data seem to favour the latter. 



1. Introduction 

A lattice investigation of the large AT limit of 
SU(TV) gauge theories (l|] is interesting in many 
respects: it can improve our understanding of 
SU(3) gauge theory as a base for full QCD and 
allows a direct comparison with predictions de- 
rived from theories beyond the Standard Model, 
often rigorous only when TV — > oo. 

Existing lattice results suggest that SU(2) and 
SU(3) may be close to the large TV limit. For this 
reason, we have investigated SU(TV) gauge theo- 
ries in 3+1 dimensions for TV = 2, 3, 4, 5 ||. If as 
TV — > oo a smooth limit exists, corrections to this 
limit can be expressed as a power series in 1/TV 2 . 
Of course, a priori it is unknown whether observ- 
ables at small TV can be continued analytically to 
TV = oo and how many terms of the series are 
needed to relate quantities at finite TV to their 
TV = oo values. To answer those questions, we 
have studied the dependence on TV of the funda- 
mental string (Sect. ^), of the mass of the ++ 
and 2 ++ glueballs (Sectj3f) and of the topological 
susceptibility (Sect. ^) g]. 

As soon as TV = 4 other strings appear in the 
spectrum that are stable (|] . These strings are de- 
noted as fc-strings, k = 1, TV/2 being an integer 
identifying the class of the representations of the 
gauge group under which sources connected by a 
£;-string transform. The values of the ratios of 
the string tensions 07./07 (the denominator be- 



ing the tension of the fundamental string) give 
insights on the mechanism of colour confinement. 
A numerical investigation of those ratios |^|,|| is 
reported in Sect. pi 

2. The string tension 

In our calculations we have used the Wilson ac- 



tion5 = /3E M >,,i(l 

with (3 = 2N/g 2 , g being the coupling constant. 

The string tension has been extracted by look- 
ing at the decay of correlation functions of 
blocked Polyakov loops. To the lowest order in 
1/L, the string tension a is related to the small- 
est mass governing this decay by the formula 
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where L is the length of the string and D is the 
number of dimensions. If the infrared behaviour 
of the string can be described by an effective 
string theory, the coefficient c s is universal || 
and depends only on the number of bosonic and 
fcrmionic modes propagating along the string. A 
reliable extraction of the string tension must keep 
into account the term proportional to 1/L. For 
this reason, we have determined c s for the gauge 
group SU(2) in D = 2 + 1 and D = 3 + 1 §. Our 
fits give c s = 1.066 ± 0.036 and c s = 0.94 ± 0.04 
respectively. Those values are compatible at the 
2cr level with an effective string theory of the 
bosonic type and are incompatible with other 
simple string models. In the following we will 
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Figure 1. The masses of the lowest-lying glueball 
states as a function of 1/N 2 with the best fits. 



assume c s — 1, i.e. that the long range fluctu- 
ations of the string are described by an effective 
string theory and that the string is bosonic. 

The string tension is used to set the scale of 
physical quantities at different N. Equivalent 
quantities are identified by equivalent couplings. 
In the original 't Hooft's idea, equivalent cou- 
plings in the large N limit correspond to fixed 
A = g 2 N. We have verified that the string ten- 
sions for different N lie on an universal curve as 
a function of A/ = q 2 N , g 2 being the tadpole im- 
proved coupling As a guidance, equivalent 
(3 at different N can be identified by using the 
relationship j3(N)/(3(N') = N 2 /N' 2 . 

3. The mass spectrum 

With the method described in we have de- 
termined the mass of the ++ and 2 ++ glueballs 
at various N. Results are available also for the 
++ *, but they are preliminary. 

Our analysis has firstly involved an extrapola- 
tion to the continuum limit for each TV, then an 
extrapolation to the large N limit of the contin- 
uum results. In Fig. |]we plot the masses in units 
of yfd of the above states as a function of 1/N 2 . 
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Figure 2. The ratio Xt as a function of 

1/N 2 with the best fit. 



A fit of the form m(N) = M + a/N 2 works for 
all these states and all the way down to SU(2). 
Hence for these quantities a smooth N — > oo limit 
exists and the first expected correction explains 
the finite N results all the way down to SU(2). 
The N — oo values extracted for M = M/y/a 
areQ M = 3.380(69) for the 0++, M = 4.94(13) 
for the 2++ and M = 6.74(16) for the 0++*. 

4. Topological susceptibility 

Large N arguments were the framework used to 
derive the Witten-Veneziano formula [pj, which 
relates the topological susceptibility Xt to the 
mass of the rjf. This formula does work even 
when the large limit of the susceptibility is 
replaced by its SU(3) value. This fact possibly 
implies that SU(3) is close to SU(oo). In Fig. || 
we report results for x¥^ IV^i which has discreti- 
sation artifacts proportional to a 2 a. Our fit gives 
xl /4 /Vv = 0.3739(59) + 0.439/A f2 for 2 < < 5. 

As A'' increases and for large values of the lat- 
tice size large correlations appear in the time se- 

1 In this work, our error bars refer to a confidence level of 
68%, while in Ref. fcf a confidence level of 75% was taken. 
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Table 1 

<Tk/crf in SU(iV) in D dimensions at zero temper- 
ature, for the values of D, N and k shown. 



D 


N 


k 




3+1 


4 


2 


1.357(29) 


3+1 


5 


2 


1.583(74) 


2+1 


4 


2 


1.3548(64) 


2+1 


6 


2 


1.6160(86) 


2+1 


6 


3 


1.808(25) 



ries of the topological charge. This phenomenon 
can be explained by the suppression of instantons 
of small size at large N, whose distribution as a 
function of N and of the size p is expected to 
behave as p ( - 11N / 3 )~ 5 Our data confirm this 
behaviour. 

5. fc-strings 

We have studied the ratio Vk/vf in both 
D=3+l and D=2+l Our results for its 

zero temperature value at various N and k are 
reported in Table 1. Those values should be com- 
pared to predictions based on various scenarios 
for the confining vacuum. Our data immediately 
exclude unbound strings (cr^ oc fee/); also the 
bag model g (cr fc cx y/[k(N - K)(N + 1)]/2N) 
is excluded. Our results in D=3+l are in agree- 
ment at the 2a level with both the Casimir scal- 
ing g (ok/a f = [k(N - k)]/(N - 1)) and the 

MQCD §] (ofc/07 = sin(fc7r/A0/sin(7r/iV)) Pre- 
dictions. Other calculations for D=3+l SU(6) 
find a better agreement with MQCD [M. Our 
results in D=2+l seem to favour Casimir scaling. 

At high temperature the spatial fc-string ten- 
sions are found to be closer to Casimir scaling 
than to MQCD (the evidence being more strik- 
ing in D=2+l). Our D=3+l high T results for 
the spatial /c-string tensions agree with a re- 
cent model calculation Jll[]. Other evidences for 
Casimir scaling in D=2+l are discussed in ||. 

6. Conclusions 

We have reported results from the first modern 
lattice calculation aimed at investigating the large 



N limit of SU(iV) gauge theories in D=3+l. We 
have shown that this limit seems to exist. More- 
over, the values of observables like the masses of 
the lowest-lying glueballs and the topological sus- 
ceptibility at finite N are related to their N = 00 
counterpart by a simple 1/N 2 correction. 

We have also investigated the ratios <Jk/&f for 
stable fc-strings. While those ratios seems to fulfill 
Casimir scaling in D=2+l in various regimes and 
possibly in D=3+l for the high T spatial string 
tensions, our data can not resolve between that 
prediction and the one based on MQCD at zero 
temperature in D=3+l. 
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